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Abstract 

In recent years answer set programming has been extended to deal with multi-valued 
predicates. The resulting formalisms allows for the modeling ol continuous problems as 
elegantly as ASP allows for the modeling of discrete problems, by combining the stable 
model semantics underlying ASP with fuzzy logics. However, contrary to the case of clas- 
sical ASP where many efficient solvers have been constructed, to date there is no efficient 
fuzzy answer set programming solver. A well-known technique for classical ASP consists of 
translating an ASP program P to a propositional theory whose models exactly correspond 
to the answer sets of P. In this paper, we show how this idea can be extended to fuzzy 
ASP, paving the way to implement efficient fuzzy ASP solvers that can take advantage of 
existing fuzzy logic reasoners. To appear in Theory and Practice of Logic Programming 
(TPLP). 
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1 Introduction 



Answer Set Programming (ASP), see e.g. (Baral 2003) is a form of non-monotonic 
reasoning based on the stable model semantics for logic programming ( Gelfond and 
Lifschitz 1988). Intuitively, in answer set programming one writes a set of rules 
(the program) such that certain minimal models (the answer sets) of this program 
correspond to solutions of the problem of interest. 

In recent work, logic programming has been extended to handle many different 



facets of imperfect information. Most notably are the probabilistic ( Baral et al. 2007 



Damasio and Pereira 2000| Fuhr 2000[ Lukasiewicz 1998||Lukasiewicz 1999[|Ng and 


Subrahmanian 1993| |Ng and Subrahmanian 1994 


Straccia 2008 


) and possibilistic 


(Alsinet et al. 2002| |Bauters et al. 2010| |Nicolas et al. 2005 


Nicolas et al. 2006 1 ex- 


tensions to handle uncertainty, the fuzzy extensions (|Cao 2000 ; 


Ishizuka and Kanai 


1985 Lukasiewicz 2006; Lukasiewicz and Straccia 2007a; Lukasiewicz and Straccia 


2007b Madrid and Ojeda-Acicgo 2008||Madrid and Ojcda-Acicg 


o2009 Saad 2009a 


Straccia 2008| Van Nieuwenborgh et al. 2007b Vojtas 2001| Wagner 1998 ) which 



allow to encode the intensity to which the predicates are satisfied, and, more gen- 
erally, many- valued extensions ( |Damasio et al. 2004| |Damasio et al. 2007| |Damasio| 
and Pereira 2001a Damasio and Pereira 2001b; Damasio a nd Pereira 2004[ |Emden| 



1986 Fitting 1991 Kifer and Li 1988 Kifer and Subrahmanian 1992[ Lakshmanan 



1994| [Lakshma nan and Sadri 1994] |Lakshmanan and Sadri 1997 Lakshman an and 
Shiri 2001 Lakshmana n 1997] |Loyer and Straccia 2002| |Loyer and Straccia 2003 



Nerode et al. 1997| [Shapiro 1983| [Straccia 2005[ |Straccia 2006] |Straccia et al. 2009 



Subrahmanian 1994 1. In this paper we focus on a fuzzy extension of ASP, called 
fuzzy answer set programming (FASP), which combines the stable model semantics 
for logic programming with fuzzy logics. More generally, FASP provides a seman- 
tics for logic programs in which the truth of predicates (or propositions) may be 
graded. Such grades may mean different things in different applications, but often 
they are related to the intensity to which a given property is satisfied. From an 
application point of view, this is important because it allows to describe continuous 
phenomena in a logical setting. Thus a formalism is obtained in which problems 
with continuous domains can be modeled with the same ease by which discrete 
problems can be modeled in classical ASP. 

In recent years, efficient solvers for classical ASP have been developed. Some of 



these are based on the DPLL algorithm (Davis and Putnam 1960) such as Smodels 



(Simons 2000) and DLV (Leone et al. 2006), others use ideas from SAT solving 



such as clasp ( Gebser et al. 2009 ) , while still others directly use SAT solvers to find 



answer sets, e.g. ASSAT (Lin and Zhao 2004), cmodels (Giunchiglia et al. 2004) 



and pbmodels (Liu and Truszczyhski 2005). The SAT based approaches have been 



shown to be fast, and have the advantage that they can use the high number of 
efficient SAT solvers that have been released in recent years. The DPLL based 
solvers have the advantage that they allow a flexible modeling language, since they 
are not restricted to what can directly and efficiently be translated to SAT, and 
that they can be optimized for specific types of programs. 



Probabilistic ASP can be reduced to classical SAT (Saad 2009b), allowing imple 
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mentations using regular SAT solvers. Likewise, possibilistic ASP can be reduced 



to classical ASP (Nicolas et al. 2006), which means ASP solvers can be used for 



solving possibilistic ASP programs. 

In the case of fuzzy ASP programs with a finite number of truth values, it has 



been shown in (Van Nieuwenborgh et al. 2007a) that FASP can be solved using reg- 



ular ASP solvers. Unfortunately, to date, no fuzzy ASP solvers or solving methods 
have been constructed for programs with infinitely many truth values. Our goal in 
this paper is to take a first step towards creating such efficient solvers by showing 
how the idea of translating ASP programs to SAT instances can be generalized 
to fuzzy answer set programs. In this way we can create fuzzy answer set solvers 
that use existing techniques for solving fuzzy satisfiability problems, e.g. based on 
mixed-integer programming or other forms of mathematical programming. Specifi- 



cally, we focus on the ASS AT approach introduced in (Lin and Zhao 2004). While 



translating ASP to SAT is straightforward when ASP programs do not contain 
cyclic dependencies, called loops, careful attention is needed to correctly cover the 
important case of programs with loops. The solution presented by ASSAT is based 
on constructing particular propositional formulas for any loop in the program. In 
this paper, we pursue a similar strategy where fuzzy loop formulas are used to 
correctly deal with loops. Our main contributions can be summarized as follows: 



1. We define the completion of a fuzzy answer set program in the sense of (Van 



Nieuwenborgh et al. 2007b I, and show that the answer sets of FASP programs 



without loops are exactly the models of its completion. 



2. By generalizing the loop formulas from (Lin and Zhao 2004), we then show 



how the answer sets of arbitrary FASP programs can be found, provided 
that the fuzzy logical connectives are t-norms. We furthermore show how 
the ASSAT procedure, which attempts to overcome the problem with an 
exponential number of loops, can be generalized to the fuzzy case. 



We furthermore show that the FASP semantics in terms of unfounded sets (Van 



Nieuwenborgh et al. 2007a ) coincide with the FASP semantics in terms of fixpoints 



(see e.g. (Lukasiewicz 2006)). This is necessary because the development of loop 



formulas can more easily be done using the unfounded semantics, while the gener- 
alization of the ASSAT procedure is based on the fixpoint semantics. 

The structure of the paper is as follows. Section [2] recalls the basic fuzzy logic 
operators and Section|3]recalls the FASP framework from (Van Nieuwenborgh et al. 
2007b ). In Section[3]we furthermore show that the unfounded semantics and fixpoint 
semantics of FASP coincide. Next, we define the completion of a FASP program in 
Section[I]and discuss the problems that occur in programs with loops. Section[5]then 
shows how these problems can be solved by adding loop formulas to the completion. 
We illustrate our approach on the problem of placing a set of ATM machines on 
the roads connecting a set of cities such that each city has an ATM machine nearby 
in Section [6} The reason for restricting our approach to t-norms is discussed in 
Section [7j Afterwards, in Section [5J we give an overview of related work and then 
present the conclusions in Section [9] 



A preliminary version of this paper appeared in ( Janssen et al. 2008). This paper 
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extends our earlier work by adding proofs, a detailed use case and a discussion 
on the problems that occur when programs are allowed to contain t-conorms in 
the body. Furthermore we improved the presentation by removing the aggregation- 
based approach that was used in the aforementioned work. 



2 Preliminaries 

In general, fuzzy logics are logics whose semantics are defined in terms of variables 
that can take a truth value from the unit interval [0, 1] instead of only the values 
(false) and 1 (true). Different ways exist to extend the classical logic connectives, 



leading to different logics with different tautologies and axiomatizations (Hajek 



2001 Novak et al. 1999). We briefly recall the most important concepts related to 
fuzzy logic connectives. 

A negator is a decreasing [0, 1] — > [0, 1] mapping A/" satisfying A/"(0) = 1 and 
Af(l) — 0. A negator is called involutive iff Vx G [0, 1] : Af(Af(x)) = x. 

A triangular norm (t-norm) is an increasing, commutative and associative 
[0, l] 2 — > [0, 1] operator T satisfying Vx € [0, 1] : T(l,x) = x. Intuitively, this op- 
erator corresponds to logical conjunction. In this paper, we restrict ourselves to 
left-continuous t-norms. As the most commonly used t-norms obey this restriction, 
this poses no practical constraint. 

A triangular conorm (t-conorm) is an increasing, commutative and associative 
[0, l] 2 — > [0, 1] operator S satisfying Vx € [0, 1] : 5(0, x) = x. Intuitively, it corre- 
sponds to logical disjunction. 

An implicator X is a [0, l] 2 —> [0, 1] operator that is decreasing in its first and 
increasing in its second argument, and satisfies 1(0,0) = 1(0,1) = 1(1,1) = 1, 
X(1,0) = and Vx £ [0, 1] : X(x, 1) = x. Every left-continuous t-norm induces a 
residual implicator defined by I(x,y) — sup{A <= [0, 1] | T(x, A) < y}. Further- 
more, a left-continuous t-norm T and its residual implicator X satisfy the residu- 
ation principle, i.e. for x,y,z in [0,1] we have that T(x,y) < z iff x < X(y,z). 
For any left-continuous t-norm T, its residual implicator X satisfies 

X(x,y) = liSx<y (1) 

For a given implicator I its induced negator is the operator Af defined by Af(x) = 
I(x, 0). We summarize some common t-norms, t-conorms, residual implicators, and 
induced negators in Tables [I] and [2] 

The biresiduum of a residual implicator X is denoted as ~, and defined by 
x sa y — T(X(x,y),I(y,x)). Note that due to ([I]) it follows that x w y is always 
equal to either X(x,y) or X(y,x). We denote the particular choice of t-norm and 
implicator using a subscript, when it is not clear from the context, as in x sa m y = 
T m {X m {x,y),X m {y,x)). 

A fuzzy set A in a universe X is an X —> [0,1] mapping. For x € X we 
call A(x) the membership degree of x in A. For convenience we denote with 
A = {a^ 1 , . . . , a^" } the fuzzy set A satisfying A(etj) = fcj for 1 < i < n and A(a) = 
for a {ai, . . . , a n }. We use ^(X) to denote the universe of all fuzzy sets in X. The 
support of a fuzzy set A is defined by supp(A) = {x e X \ A(x) > 0}. Inclusion 
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t-norm t-conorm 



T m (x,y) = mm(x,y) 

Ti(x, y) = max(0, x + y — 1) 

T P {x,y) = x-y 



S m (x,y) = max(i,i/) 
Si(x,y) = mm(x + y, 1) 
S P (x,y) = x + y~x-y 



Table 1. Common fuzzy t-norms and t-conorms over ([0, 1], <) 



t-norm residual implicator 



induced negator 



Tin 

77 



Xm{x,y) 



T p {x,y) 



y if x > y 
1 otherwise 
Xi (x, y) = min(l, 1 — x + y) 
y/x if x > y 
1 otherwise 



A/" m (z) 



if x > 

1 otherwise 
JV; (a;) = 1 - x 



N P {x) =N m (x) 



Table 2. Common residual pairs and induced negators over ([0, 1], <) 



of fuzzy sets in the sense of Zadeh is defined as A C B iff Vx € X : < B(x). 

Last, in this paper we will write the difference A Q B of two fuzzy sets to denote 
the fuzzy set defined by (A Q B)(x) = max(0, A(x) — B(x)). 

A signature is a tuple § = (A, T, C, I, N) , with A a set of atoms (i.e. propositional 
letters), T a set of t-norms, C a set of t-conorms, I a set of implicators, and N a 
set of negators. Additionally, we demand that max and min are in the signature. A 
fuzzy formula over a signature S then is either an atom, a value from [0, 1], or the 
application of a t-norm or t-conorm from T, resp. C, to two formulas, the application 
of an implicator from I to two formulas, or the application of a negator from N to a 
single formula. A fuzzy theory over a signature S is a set of fuzzy formulas over S. 
An interpretation I over a signature S is a mapping from A to [0, 1]. It is extended 
to fuzzy formulas in a straightforward way, i.e. if F and G are fuzzy formulas, 
then I(T(F,G)) = T(I(F),I(G)) (with T € T), I(S(F,G) = S (1(F), 1(G)) (with 
S e C), I(2(F, G)) = 1(I(F), 1(G)) (with 1 e I), and I(Af(F)) = Af(I(F)) (with 
Af € N). An interpretation M is a model of a fuzzy formula F, denoted M \= F, iff 
M(F) = 1. An interpretation M is a model of a fuzzy theory 0, denoted M |= 0, 
iff for each F € 9 we have that M |= F. 

A particular signature leads to a particular fuzzy logic (Hajek 2001). For ex- 
ample, S = (A, {7m}, {S m }, {l m }, Wm}) gives rise to Godel logic, S = (A, 
{77,7^}, {Si,S m }, {%}, {Mi}) gives rise to Lukasiewicz logic, and§ = (A, {7^,7™}, 
{S p ,S m }, {Ip}, {Np}) gives rise to product logic. For example, Lukasiewicz logic 
is generally considered to be closest in spirit to classical logic, in the sense that 
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many of its important properties are preserved. Another important advantage of 
Lukasiewicz logic is that the implicator is continuous, which is not the case for Godel 
or product logic. Reasoning in this logic can be done using mixed integer program- 
ming, whereas reasoning in Godel logic can be done with the help of boolean SAT 
solvers. As in the boolean case, satisfiability checking in these three particular logics 
is NP-complete. 

Last, we denote the infimum, resp. supremum of two elements of [0, 1] as a n b, 
resp. aUb. 



3 Fuzzy Answer Set Programming 

Over the years many different fuzzy answer set programming formalisms have 



been developed (Lukasiewicz 2006 Lukasiewicz and Straccia 2007a Madrid and 



Ojeda-Aciego 2008[ |Madrid and Ojeda-Aciego 2009[ |Saad 2009a] |Straccia 2008} 
Van Nieuwenborgh et al. 2007b). Most of these base their semantics on fixpoints 



or minimal models, in combination with a reduct operation. The approaches de- 



scribed in (Loyer and Straccia 2006 Van Nieuwenborgh et al. 2007b I, however, are 



constructed from a generalization of unfounded sets. As the development of loop 
formulas can be done more elegantly when starting from unfounded sets, and we 
can show that the fixpoint semantics are equivalent to the semantics proposed in 



(Van Nieuwenborgh et al. 2007b), we will base our semantics on the latter frame- 



work. However, because the generalization of the ASS AT procedure from (Lin and 



Zhao 2004) is based on fixpoint semantics, in this section we also show novel links 



between the unfounded and fixpoint semantics that ensure the correctness of our 



generalized procedure. First, we recall the main definitions from (Van Nieuwenborgh 



et al. 2007b[) 

A litera^J is either an atom a or a constant from [0, 1]. An extended literal is 
either a literal (called a positive extended literal) or an expression of the form 
J\f(a) (called a negative extended literal) , with a an atom and Af an arbitrary 
negator. A rule r is of the form 

r : a T(b 1: ...,b n ) 

where n > 0, a is a literal, {b%, . . . , b n } is a set of extended literals, T is an arbitrary 
t-norm, and r is a rule label. Furthermore, for convenience, we define Tib) = b, and 
define 7~(&i, ■ ■ ■ , b n ) as T(b\, T{b%, ■ ■ ■ ))■ The literal a is called the head, H(r), of 
r, while the set {&i, . . . , b n } is called the body, B(r), of r. We use Lit + (B(r)) to 
denote the set of positive extended literals from the set B{r). Given a rule r we 
denote the t-norm used in its body as T r ] the residual implicator corresponding to 
T r is denoted as I r . A constraint is a rule with a constant in its head, whereas 
a fact is a rule with a constant as its body (and no constant in its head). For 
convenience, we abbreviate a rule of the form r : a <— T(b, 1), with b an extended 
literal, as r : a <— b. 

A FASP program P is a finite set of rules. We call a program simple if no rule 



1 As usual, we assume that programs have already been grounded. 
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contains negative extended literals. The set of all atoms occurring in P is called the 
Herbrand Base Bp of P. Note that the Herbrand Base is finite since we assume 
that no function symbols occur in the terms of the ungrounded program. For any 
a £ Bp we define the set P a as the set of rules with atom a in their head. An 
interpretation / of P is a fuzzy set in Bp. We extend interpretations to constants 
from [0, 1], extended literals, and rules as follows: 

1. 1(c) = c if c £ [0,1] 

2. I(Af(l)) = Af(I(l)) if / is a literal 

3. I(a^T{bi,...,b ri ))=X r (T{I{b 1 ),...,I{b ri )),I(a)) 

A model of a program P is an interpretation / of P such that for each rule r £ P 
we have I(r) = 1. 

Note that, although each rule in a FASP program can only have a single t-norm 
in its body, a rule with mixed t-norms, such as r : a <— 7~i(a, T^ib, c)), can easily 
be simulated by introducing a polynomial number of new literals and rules. In the 
case of rule r we need one new literal a' and two new rules r% : a <— 7i(a, a!) and 
r 2 : a' <- Ti(b, c). 

Example 1 

Consider the program P, which consists of the following set of rules: 

rf.ai- T m (b,c) 
r 2 :b^- 0.8 

r 3 :ci-T m (a,M(b)) 
r 4 : <- Ti(a, b) 

Note that rule r 2 is a fact, and rule r 4 a constraint. The fuzzy sets Ji = {a , 6 8 , c }, 
and I2 — {a 2 , b S ,c 2 } are interpretations of P. For both interpretations we have 
that h(ri) — hfri) = . . . = hfa) = h{ri) = 1, i.e. they both are models of the 
program. 

The definition of fuzzy answer sets relies on the notion of unfounded sets, stud- 



ied in (Van Nieuwenborgh et al. 2007b I, which correspond to sets of "assumption 



atoms" that have no proper motivation from the program. 



Definition 1 (Unfounded sets (Van Nieuwenborgh et al. 2007b]) ) 



Let P be a FASP program and let / be an interpretation of P. A set U C Bp is 
called unfounded w.r.t. I iff for each atom u £ U and rule r £ P u either: 

1. Ur\Lit+{B(r)) ^ 0; or 

2. I(u) > I(B(r)); or 

3. I(B(r)) = 

An interpretation I of P is called unfounded-free iff supp(I) f)U = for any set 
U that is unfounded w.r.t. /. 



Intuitively, an unfounded set w.r.t. an interpretation / of a FASP program P is 
a set of atoms that obtain a value in / that is not motivated by the rules of the 
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program. The first condition of Definition [l] ensures that the values of the literals 
in U are justified by the values of literals not in U. The second condition shows 
that the degree to which a rule can motivate an atom is bounded by the value of 
its body. The third condition is needed to obtain a proper generalization of the 



classical definition of unfounded sets (Baral 20031 (see (Van Nieuwenborgh et al. 



2007b ) for more details) . 
Example 2 

Consider program P and interpretations I\ and I 2 from Example [I] For I 2 we can 
see that U 2 = {ct, c} is an unfounded set, as for rule r\ and r^, the only rules with 
a or c in the head, we have that Lit + (B(r 1 )) C\U 2 ^$ and Lit + (B(r 3 )) C\U 2 ^$. 
Since suppil^) n U 2 ^ 0, interpretation I 2 is not unfounded-free. Interpretation Ii, 
however, is unfounded-free. 



As shown in (Van Nieuwenborgh et al. 2007b), answer sets of FASP programs 



can be defined as the unfounded-free interpretations, which reflects the intuition 
that each atom in an answer set should have a proper motivation from the program. 



Definition 2 [Answer Set ( Van Nieuwenborgh et al. 2007b)) 



Let P be a FASP program. A model M of P is called an answer set iff M is 
unfounded-free. 

Example 3 

Consider program P and interpretation 1\ from Example [I] Since we know from 
Example [2] that I\ is unfounded-free, it follows that I\ is an answer set of P. 

An alternative definition of answer sets, in terms of fixpoints, exists (see for 



example (Lukasiewicz 2006)). We will use this to generalize the ASSAT procedure 
described in (Lin and Zhao 2004). 



Definition 3 (Immediate Consequence Operator (Lukasiewicz 2006)) 



Let P be a FASP program. The immediate consequence operator of P is the 

mapping Up : J- {Bp) —> J- (Bp) defined by 

U P (I)(l) = sup{I(B(r)) \ rePi} 



As shown in (Lukasiewicz 2006), for simple programs this operator is monotonic 



and thus has a least fixpoint (Tarski 1955), denoted as rip. For these simple pro 



grams, (Lukasiewicz 20061 then defines the answer sets of a program as the least 



fixpoints of this operator. Since this operator is monotonic, the least fixpoint is 
unique and can be found by iteratively applying lip from the interpretation until 



a fixpoint is encountered. For non-simple programs, (Lukasiewicz 2006) defines a 



reduct operation that transforms a non-simple program into a simple program. 
Definition 4 (Reduct (Lukasiewicz 200~6 )) 



Let P be a FASP program and let r : a <— T(b\, . . . , b m ,J\f(b m+1 ), . . . ,J\f(b n )) be a 
rule in P, where (b\, . . . , b m ) = Lit + (B(r)). The reduct of rule r, with respect to 
an interpretation /, is denoted as r 1 , and defined by 

r 1 T(h, b m , I(Af(b m+1 )), I(Af(b n ))) 

The reduct of a program P w.r.t. an interpretation / is defined as P 1 — {r 1 | r G P}. 
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Example 4 

Consider program P and interpretation I\ from Example [l] The reduct of P with 
respect to 1\ then is the following program 

rl 1 : a <- T m (b, c) 
r^ 1 : 6 <- 0.8 
r^ 1 : c^7;„(a,0.2) 
r^ 1 : <- 7I(a, 6) 

In the following we show the novel result that the semantics in terms of fixpoints 
coincide with those in terms of unfounded sets. An important lemma regarding the 
immediate consequence operator and reduct is the following. 

Lemma 1 

Let P be a FASP program. For any interpretation / of P it holds that L = Hp(I) 
xSI = IL P i(I). 

Proof 

Follows trivially by the construction of P 1 and Definition |ZJ □ 

We now show that any answer set is a fixpoint of the immediate consequence 
operator. 

Lemma 2 

Let P be a FASP program. Then any answer set A of P is a fixpoint of lip. 
Proof 

Let A be an answer set of P. We show that A(a) — sup{A(B(r)) \ r 6 P a } = 
Hp(A)(a) for any a £ Bp, from which the stated readily follows. The proof is split 
into the case for a G supp(A) and a (jL supp(A). For any a € supp(A) it must 
hold that {a} is not unfounded w.r.t. A, meaning that P a ^ and there is some 
r e P a such that A(a) < A(B(r)). Since A(r) = 1, it then follows from that 
A(a) = A{B{r)). As for any r' € P a we have A(r') — 1, from ([lj it also follows 
that A(B(r)) = A(a) > A(B(r')). Hence A{B(r)) = A(a) is the supremum of 
{A(B{r')) | r' € P a }. 

The case for a $ supp(A) is as follows. First remark that as A(a) = 0, it follows 
from and the fact that A{r) = 1 for each r € P a , that A(B(r)) = 0. Hence, 
A(a) = sup{A(B{r)) \ r € P Q }. □ 

Second we show that answer sets can be characterized in terms of fixpoints of 
the immediate consequence operator. 

Proposition 1 

Let P be a FASP program. An interpretation A is an answer set of P iff A = H* pA . 
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Proof 

Let M be a model of P. In (Van Nieuwenborgh et al. 2007b) it was shown that the 



least fixpoint of H p m must necessarily be unfounded-free (Proposition 4). As any 
fixpoint of lip is a model of P, we only need to show that if M is unfounded- free, 
it is the least fixpoint of UpM . 

Suppose M ^ Tl*p M ■ Then, since any unfounded-free model is a fixpoint of IIpM 
due to Lemmas[l]and [2j it holds that some set M' C M exists such that M' = U* pM ■ 

Consider then U = {u £ B P \ M'(u) < M(u)}. Surely U C supp(M) and hence 
U H supp(M) 7^ 0. We now show that U is unfounded with respect to M, leading 
to a contradiction. First, we show that for any atom u £ U and rule r £ P u it holds 
that 

(lAt + (B(r)) nU = |)W M(B(r)) < M{u) (2) 

as follows 

Lit+(B(r)) n U = 
= (Def. n) fil £ Lit+(B(r)) : I £ U 

= (Duality V,3> £ Lit+(B(r)) : I & U 

= (Def. U) VI £ Lit+(B(r)) : M(a) = M'(a) 

=> (M(l) = M(l M )) M(B(r M )) = M'{B{r M )) 

(Monotonicity sup) M{B(r M )) < swp r , eP M M'(B(r')) 
= (M' = TlpM (M')) M(B{r M )) < M'(u) " 
=> (u G U, Def. U) M{B(r M )) < M(u) 

Thus, since it follows from the Definition of r M that M(B(r M )) = M(B(r)), we 
have shown that ^ holds. From this equation we obtain that 

(Lit + (B{r)) n U ^ 0) V (M(B(r)) < M(ufj 

Hence 

(Lit + {B(r)) n U ^ 0) V (M(B(r)) < M(uj) V (M(B(r)) = Oj 
Which means U is unfounded with respect to M, a contradiction. □ 
Example 5 

Consider program P and interpretations I\ and I2 from Example [T] Computing the 
least fixpoint of ilp/j and IIpj 2 can be done by iteratively applying these operators, 
starting from 0, until we find a fixpoint. Hence for P Jl we obtain in the first iteration 
Ji = n P / 1 (0) = {b - 8 }. The second iteration gives J 2 = 11^ (Ji) = {b 8 } = J u 
hence a fixpoint, meaning {b o s } = 1\ is the least fixpoint of Hp^. Iteratively 
applying n p j 2 brings us J\ = IIp/ 2 (0) = {b 8 }, which is also a fixpoint of IIp/ 2 . 



4 Completion of FASP programs 

In this section we show how certain fuzzy answer set programs can be translated 
to fuzzy theories such that the models of these theories correspond to answer sets 
of the program and vice versa. Such a correspondence is important as it allows us 
to find answer sets using fuzzy SAT solvers. 
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Definition 5 [Completion of a FASP program) 

Let P be a FASP program. The completion of P, denoted as comp(P), is defined 
as the following set of fuzzy formulas: 

{a ps (max{S(r) \ r £ P a }) \ a £ B P } U {l r (B(r),H(r)) \ r £ P, H(r) € [0, 1]} 

where ss is the biresiduum of an arbitrary residual implicator, and T r is the residual 
implicator of the t-norm used in the body of rule r. 

The completion of a program consists of two parts, viz. a part for the literals 
{a ss (max{£?(r) | r £ P a }) | a £ Bp}, and a part for constraints {X r (B(r), H(r)) \ 
r £ P,H(r) £ [0,1]}. The constraints part simply ensures that all constraints 
are satisfied. The literal part ensures two things. By definition of the biresiduum 
and the fact that I(a, b) = 1 iff 1(a) < 1(b) for any residual implicator, we have 
that I (a ps b) = 1 iff 1(a) < 1(b) and 1(b) < 1(a). Hence, the literal part of the 
completion establishes that rules are satisfied and second that the value of the literal 
is not higher than what is supported by the rule bodies. 

Example 6 

Consider program P from Example [T] Its completion is the following set of fuzzy 
propositions 

a ps T m (b, c) 
b ps 0.8) 

c Ps T m (a,Af m (b)) 
Xi{Ti{a,b),Q) 

Note that when applying Definition [5] for a literal I that does not appear in the 
head of any rule, we get a ps max0, where we define max0 = 0. 

We can now show that any answer set of a program P is a model of its completion 
comp(P). 

Proposition 2 

Let P be a FASP program and let comp(P) be its completion. Then any answer 
set of P is a model of comp(P). 

Proof 

Suppose A is an answer set of P. By Lemma[2] it follows that A is a fixpoint of lip, 
hence for each a £ Bp, A(a) — svp{A(B(r)) | r £ P a }- By construction of comp(P) 
and the fact that A is a model of P, it then easily follows that A |= comp(P). □ 

Example 7 

Consider program P and interpretation I\ from Example [l] It is easy to see that 
1\ is a model of comp(P). 

The reverse of Proposition [2] is not true in general, which is unsurprising because 
it is already invalid for classical answer set programming. The problem occurs for 
programs with "loops" , as shown in the following example. 
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Example 8 

Consider program P and interpretation I2 from Example [T] We can easily see that 
I2 is a model of comp(P), but, as we have seen in Example [2j it is not an answer 
set of P. 

One might wonder whether taking the minimal models of the completion would 
solve the above problem. The following example shows that the answer is negative. 

Example 9 

Consider the following program P m i n 

a <— a 

The completion comp(P m i n ) is 

a fa a 

Consider now the interpretation / = {a 2 ,p QA }. Since 1(a) = 1(a) and 

Ti(M(I(p)),M(I(a))) = max(0, 1 - I(p) + 1 - 1(a) - 1) = 0.4 

we can see that / is a model of comp(P m i n ). We show that it is a minimal model 
as follows. Suppose some I' C I exists. Then we can consider three cases: (i) 
I'(a) < 1(a) and V (p) = I(p); (ii) I'(a) = 1(a) and V (p) < I(p); (hi) I'(a) < 1(a) 
and I'(p) < I(p). In all three cases we obtain that Ti(M(I'(p)),M(I'(a))) > 0.4 > 
I'(p), since Mi(I'(a)) = 1 - I'(a) > 0.8 or Ni(I'(p)) > 0.6. Hence V is not a model 
of comp(P m i n ) and / is thus a minimal model of comp(P m i n ). However, I' is not 
an answer set of P m in since _{a° ,p ' 4 }. 

As in the crisp case however, when a program has no loops in its positive de- 
pendency graph, the models of the completion and the answer sets coincide. First 
we define exactly what a loop of a FASP program is, and then we show that this 
property indeed still holds for FASP. 

Definition 6 (Loop) 

Let P be a FASP program. The positive dependency graph of P is a directed 
graph G P = (B P ,D) where (a, 6) e D iff 3?- e P a : b e Lit + (B(r)). For ease of 
notation we also denote this relation with (a, b) € Gp for atoms a and b in the 
Herbrand base of P. We call a non-empty set L C Bp a loop of P iff for all literals 
a and b in L there is a path (with length > 0) from a to b in Gp such that all 
vertices of this path are elements of L. 

Example 10 

Consider program P m in from Example [9j The dependency graph of P m in is pictured 
in Figure [TJ We can see that {a} is a loop. If this loop was not in the program, its 
completion would become 

a pa 

P~T m (N m (p),N m (a)) 
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Fig. 1. Dependency graph of program P m i n from Example [9] 




Fig. 2. Dependency graph of program P from Example [T] 

This fuzzy theory has no models. Since program P m i n has no answer sets, this 
means the answer sets coincide with the models of the completion when removing 
the loop. 

Example 11 

Consider program P from Example [I] The dependency graph of P is pictured in 
Figure [2j We can clearly see that there is a loop between nodes a and c. Due to 
this loop, the values of a and c are not sufficiently constrained in the completion. 

From the preceding examples one might think that removing the loops from the 
program would be sufficient to make the models of the completion and the answer 
sets coincide. However, this is not the the semantics of the program then 

changes, as illustrated in the following example. 

Example 12 

Consider program P c hange consisting of the following rules 

n : a <- 0.3 
ri : a <— b 
r 3 : b <— a 

Its single answer set is {a ' 3 , b 03 }. If we remove rule r^ or r 3 , the answer set of the 
resulting program is {a 3 }. 

We can now show that for programs without loops the answer sets coincide with 
the models of their completion. We first introduce two lemmas. 

Lemma 3 

Let G = (V, E) be a directed graph with a finite set of vertices and X C V with 
X ^ 0. If every node in X has at least one outgoing edge to another node in X, 
there must be a loop in X . 
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Proof 

From the assumptions it holds that each x € X has an outgoing edge to another 
node in X. This means that there is an infinite sequence of nodes x\, x%, . . . such 
that (xi,Xi+i) € E for i > 1. Since X is finite, it follows that some vertex occurs 
twice in this sequence, and hence that there is a loop in X. □ 

Lemma 4 

Let P be a FASP program, / an interpretation of P and U C Bp. Then if 7 |= 
comp(P) and [/ is unfounded w.r.t. J it holds that for each u in U n supp(I) there 
is some r in P u such that Lit + (B(r)) DUD supp(I) ^ 0. 

Proof 

Assume that u € Uf]supp(I), in other words u G U and J(u) > 0. As {/ is unfounded 
w.r.t. /, for each r 6 P u it holds that either 

1. Lit+(B(r)) n £7 7^ 0; or 

2. /(B(r)) < J(u); or 

3. I(B(rJ) = 

We can now show that there is at least one rule r £ P u that violates the second 
and third of these conditions, meaning it must satisfy the first. 

From / |= comp(P) we know by construction of comp(P) that for each u € U, 
I(u) = sup{J(S(r)) | r <= P u }. Hence for each u € U there is a rule r £ P u such 
that I(u) = I(B(r)), thus the second condition is violated. Since I(u) > 0, it then 
also follows that the third condition is violated. 

In other words there must be some r € P u such that I(B(r)) = I(u) and 
I(B(r)) 0. Since U is unfounded w.r.t. /, this means that Lit + (B(r)) fl U ^ 
0. Since I(B(r)) ^ implies that Lit + (B(r)) C supp(I) due to the fact that 
T(0, a;) = for any t-norm T, we can conclude that there is some r € P u such that 
Lit+{B{r)) n J7 H 7^ 0. □ 

Using these lemmas we can now show that the answer sets of any program without 
loops in its dependency graph coincide with the models of its completion. This 



resembles Fages' theorem on tight programs in classical ASP (Fages 1994). 



Proposition 3 

Let P be a FASP program. If P has no loops in its positive dependency graph it 
holds that an interpretation / of P is an answer set of P iff I \= comp(P). 

Proof 

We already know from Proposition[2]that any answer set of P is necessarily a model 
of comp(P), hence we only need to show that every model of comp(P) is an answer 
set of P under the conditions of this proposition. As I \= comp(P) , it holds that / 
is a model of P. We show by contradiction that I is unfounded-free. Assume that 
there is a set U C Bp such that U is unfounded w.r.t. I and U<~) supp(I) ^ 0. From 
Lemma [4] we know that for each u € U fl supp(I) it holds that there is some rule 
r e P u such that Lit+(B(r)) n U H supp(I) ^ 0. Using the definition of G P , this 
means that for each such u there is some v! G U n supp(I) such that Gp(u,u'). 
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This however means that there is a loop in Gp by Lemma [3j contradicting the 
assumption. □ 

Hence, finding the answer sets of a program with no loops in its positive depen- 
dency graph can be done by finding models of its completion. 



5 Loop Formulas 

As mentioned in the previous section, sometimes the models of the completion are 
not answer sets. In this section, we investigate how the solution that has been 
proposed for boolean answer set programming, viz. adding loop formulas to the 



completion (Lin and Zhao 2004), can be extended to fuzzy answer set programming. 

For this extension, we start from a partition of the rules whose heads are in some 
particular loop L. Based upon this partition, for every loop L we define a formula 
in fuzzy logic, such that any model of the completion satisfying these formulas is 
an answer set. 

For any program P and loop L we consider the following partition of the rules 
in P whose head belongs to the set L (due to ( |Lin and Zhao 20"04 )) 



Rp-(L) ={s^5|((fl^B)eP)A(ffl€L)A(B+nl^)} (3) 
R, p {L) = {a <- B | ({a <- B) E P) A (a E L) A (B+ n L = 0)} (4) 

Note that this partition only takes the positive occurrences of atoms in the loop 
into account. Intuitively, the set Rp~(L) contains the rules that are "in" the loop L, 
i.e. the rules that are jointly responsible for the creation of the loop in the positive 
dependency graph, whereas the rules in Rp (L) are the rules that are outside of this 
loop. We will refer to them as "loop rules" , resp. "non-loop rules." 

Example 13 

Consider program P from Example [lj It is clear that for the loop L — {a, c} the 
set of loop rules is Rp~(L) — {r\,r$\ and the set of non-loop rules is Rp(L) = 0. 

Example 14 

Consider program P from Example [T] with interpretations I\ and I2 from Examplc[T] 
once again. It is clear that in I\ no loop rules were used to derive the values of a 
and c, whereas in I2 only loop rules are used. 

Hence there is a problem when the value of literals in a loop are only derived from 
rules in the loop. To solve this problem, we should require that at least one non-loop 
rule motivates the value of these loop literals. As illustrated in the next example, 
one non-loop rule is sufficient as the value provided by this rule can propagate 
through the loop by applying loop rules. 

Example 15 

Consider program P c hange from Example 12 again. Clearly this program has a loop 



L = {a,b} with Rp(L) = {r 2 , r 3 } and Rp(L) = {r\}. Consider then interpretations 
7i = {a 3 , b 03 } and I 2 = {a 1 , b 1 }. We can easily see that I\ is an answer set of P, 
whereas 1% is not, although they are both models of comp(P). The problem is that 
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in I2 the values of a and b are higher than what can be derived from the non-loop 
rule ri, whereas in I\ their values are exactly what can be justified from applying 
rule r\. The latter is allowed, as values are properly supported from outside the 
loop, while the former is not, as in this case the loop is "self-motivating". 

To remove the non-answer set models of the completion, we add loop formulas 
to the completion, defined as follows. 

Definition 7 {Loop Formula) 

Let P be a FASP program and L = {h, . . . , l m } a loop of P. Suppose that Rp(L) = 
{ri, . . . , r„}. Then the loop formula induced by loop L, denoted by LF(L, P), is the 
following fuzzy logic formula: 

J(max(Zi, . . . , l m ), max(S(n), . . . , B{r n )) (5) 

where X is an arbitrary residual implicator. If Rp(L) = 0, the loop formula becomes 

T(max(Zi,...,Z m ),0) 



The loop formula proposed for boolean answer set programs in (Lin and Zhao 



2004) is of the form 

-.(/\ B(n) V ... V /\ B{r n )) Hi A ... A ^l m ) (6) 

It can easily be seen that ^ is a straightforward generalisation of ^ as the 
latter is equivalent to 

(h V . . . V l m ) => (/\ B(n) V ... V /\ B(r„)) 

Note that this equivalence is preserved in Lukasiewicz logic, but not in Godel or 
product logic. Furthermore, since / (= I(max(ii, . . . , / m ), 0) only when max(/(ii), 
. . . , I(l m )) < 0, it is easy to see that in the case where Rp(L) — 0, the truth value 
of all atoms in the loop L is 0. 

Example 16 

Consider program P and interpretations I\ and I2 from Example [l] The loop for- 
mula for its loop L — {a, c} is the fuzzy formula I m (max(a, c), 0), since Rp(L) = 0. 
It is easy to see that I2 does not satisfy this formula, while interpretation 1\ does. 

Example 17 

Consider program P c hange from Example 12 The loop formula for its loop L = 



{a,b} is the propositional formula X m (max(a, b), 0.3), since Rp(L) — {ri}. Again 
we see that interpretation 1\ from Example [15] satisfies this loop formula, whereas 
interpretation I 2 from the same example does not. 

We now show that by adding loop formulas to the completion of a program, we 
get a fuzzy propositional theory that is both sound and complete with respect to 
the answer set semantics. First we show that this procedure is complete. 

Proposition 4 (Completeness) 

Let P be a FASP program, let C be the set of all loops of P, and define LF(P) = 
{LF(L, P)\Le £}. For any answer set / of P, it holds that I (= LF (P) U comp(P). 
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Proof 

Suppose 7 is an answer set of P and 7 ^ LF(P) U comp(P). Since any answer set 
is a model of comp(P) according to Proposition [2j this means that 7 \/= LF(P). 
Hence, the loop formula of some loop L in P is not fulfilled; this means: 

sup7(u) > sup I(B(r)) 

u£L reRp(L) 

Consider then the set U = {u 6 L \ I(u) > sup refl -^ L ^ I(B(r))}. We show that U 
is unfounded w.r.t. 7, i.e. we show that for each u G U and rule r G P u , at least one 
of the conditions of Definition [I] applies. 

Since P u — Rp u (L) U Rp (L), each rule r G P u must either be in Rp (L) or in 
Rp (L). We consider the following cases: 

1. If r € Rp (L) then by construction of U it holds that I(B(r)) < I(u). 

2. If r £ Rp u (L) and I(B(r)) < sup r , eR - (L) 7(P(r')), by construction of t/ we 
have that I(B(r)) < I(u). 

3. Suppose r € and 7(£(r)) > sup r , efl - {L) I{B(r')). Since T(a;,y) < 
min(x, y) for each t-norm T, we know that I(B(r)) < 1(1) for each I 6 Lit + (B(r)). 
Hence for each Z G Lit + (B(r)) we have 7(Z) > sup r , eR - ^ I(B(r')). This means 

that, since r € Pp(L) and thus Lit + (B(r)) n P 7^ 0, we know from the definition 
of £7 that Lit+(B(r)) n £7 7^ 0. 

Now remark that f7 n supp(I) 7^ as [7 C supp(I) due to 7(it) > for each 
w G P. From the above we can thus conclude that P is unfounded w.r.t. 7, and 
since U n supp(I) 7^ 0, that 7 is not unfounded-free: a contradiction. □ 

Second we show that adding the loop formulas to the completion of a program 
is a sound procedure. 

Lemma 5 

Let G = (V, E) be a directed graph and X C V, with V finite, such that each 
node of X has at least one outgoing edge to another node in X. Then there is a set 
LCI such that L is a maximal loop in X and for each I G L we have that there 
is no x G X \ L for which (I, x) G P. 

Proof 

From Lemma[3]we already know that there must be a loop in X. Hence, there must 
also be a maximal loop in X. First, remark that maximal loops must of course be 
disjoint as otherwise their union would form a bigger loop. Consider then the set 
X, which is a collection of disjoint maximal loops L and remaining nodes S (single 
nodes that are not in any loop). There is an induced graph G' of G with nodes 
S U L (i.e. each maximal loop is a single node in the induced graph) and edges E 
induced as usual (i.e. (Pi,P 2 ) G E if for some node 1% in Pi there is a node I2 in 
P2 such that (h, I2) G E and likewise for the nodes in S). Clearly, G' is acyclic as 
otherwise the nodes in G' on the cycle would create a bigger loop in X. Hence, G' 
has leafs without outgoing edges. However, a leaf cannot be in S since that would 
imply a node in X without an outgoing edge. Thus we can conclude that all leafs 
in G' are maximal loops in X. □ 
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Proposition 5 (Soundness) 

Let P be a FASP program and let LF(P) be the set of all loop formulas of P. Then 
for any interpretation 7 of P it holds that if 7 |= LF(P) U comp(P), then 7 must 
be an answer set of P. 



Proof 

Suppose I \= LF(P) U comp(P) and 7 is not an answer set of P. Since any model 
of comp(P) must be a model of P, this must mean that 7 is not unfounded-free, 
i.e. that there exists a set U C Bp such that U is unfounded w.r.t. 7. From Lemma|| 
we know that for each u £ U f~l supp(I) there must be some r £ P u such that 
Lit + (B(r)) n U n supp(I) =^ 0. Hence, by definition of Gp this means that for each 
u £ U C\ supp(I) there is some u' G U n supp(I) such that (u,u') G Gp. Using 
Lemma [5] this means that there is a set L C Pn supp(I) such that 7 is a loop in P 
and for each I G 7 there is no u e (P n supp(I)) \ L such that (Z, u) G Gp. In other 
words, for each I G 7 and rule r € p we have that 

(V n supp(J) n Lit + (B(r)) ^ij^^ln Lit + (B(r)) ^ 0) (7) 

Now, consider Z G P Since L C U P\ supp(I), we know that /(7) > 0. Hence, if 
I(B(r)) = 1(1) for some rule r e p, we know that I(B(r)) > 0. As P is unfounded 
w.r.t. P it follows from Definition [l] that L n Lit + (B(r)) ^ 0. 

Using contraposition, this means that for each I £ L and r £ P we have that 

(L n Lit+(B(r)) = 0) (l(B(r)) ? 1(1)) (8) 

By the definition of comp(P), however, we know that for each model of comp(P) 
and for each a £ Bp and r G P a we have /(a) > I(B(r)). Hence for each / G L and 
r G P from ([8]) we have that 

(L n Lit+(B(r j) = 0) (/(P(r)) < 7(1)) (9) 

Now, for each / G L and r G Kp(L) n p by definition of R~p(L) it holds that 
7 H Lit+(B(r)) = 0, meaning 7(P(r)) < 7(7). Thus, sup{7(P(r)) | r £ P P (7)} < 
sup{7(7) I I £ 7}, meaning 7 ^ LF(7, P), a contradiction. □ 

A straightforward procedure for finding answer sets would now be to extend 
the completion of a program with all possible loop formulas and let a fuzzy SAT 
solver generate models of the resulting fuzzy propositional theory. The models of 
this theory are the answer sets of the program, as ensured by Propositions [4] and 
[5] As there may be an exponential number of loops, however, this translation is 
not polynomial in general. A similar situation arises for classical ASP. The solution 



proposed in (Lin and Zhao 2004) overcomes this limitation by iteratively adding 
loop formulas. In particular, a SAT solver is first used to find a model of the 
completion of a classical ASP program. Then it is checked in polynomial time 
whether this model is an answer set. If this is not the case, a loop formula, which is 
not satisfied by the model that was found, is added to the completion. The whole 
process is then repeated until an answer set is found. We will show that a similar 
procedure can be used to find answer sets of a FASP program. 
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Starting from the fixpoint characterization of answer sets of FASP programs, we 
show that for any given model of the completion that is not an answer set, we can 
construct a loop that is violated. 

Proposition 6 

Let P be a FASP program. If an interpretation I of P is a model of comp(P) and 
I i= U* pI , then some L C supp(I W pI ) must exist such that I ty= LF(P, L). 

Proof 

Suppose I is an interpretation of P and I \= comp(P), then from the definition of 
comp(P) and Lemma[IJ we can easily see that / is a fixpoint of H P i . Since / =/= n* pI , 
some I' C I must exist such that /' = H* pI . 

Consider then the set [/ = {«£ Bp \ I(u) > I'(u)}. It holds that U = supp(lQl') 
since Pel and thus U = supp(I Q H* pI ) by definition of V . From the proof of 
Proposition [T] we then also know that for this set U the following property holds 

Vu€ U: Vr € P u : (ut + (B(r)) n U = 0) => (l(B(r)) < I(uf) (10) 

We can then show that there is a loop in U whose loop formula is violated. Since 
I = Hpi(I) we know from Lemma [l] that I = IIp(J). From the definition of lip 
this means 

Vie Bp: 1(1) = sup{ I(B(r)) \ r G Pi} 
Since the supremum is attained because P is finite we obtain 

Vie Bp: 3rePr. 1(1) = I(B(r)) 

As U C Bp this means 

Vu € U : 3r € P u : 1(1) = I(B(r)) 



Using ( 10 ) it then holds that 

VueU: 3r e P u : Lit+(B(r)) C\U 7^ 

From the definition of Gp we thus get 

Vit e U: 3u' e U: (u,v!) e Gp 

Using Lemma [5] it follows that there is a set L C U that is a loop in P such that for 
each I e L there is no I' e U \ L such that (/, I') e E. In other words, for each / e L 
there is no V e U \ L such that there is a rule r e Pi for which I' e Lit + (B(r)). 
Hence for each I e L and rule r e Pi such that U H Lit + (B(r)) ^ 0, it follows that 



L n Lit + (B(r)) 7^ 0. From (10 1 and using contraposition this means there is some 
L CU that is a loop in P and for each I e L and r e Pi if L n Lit + (B(r)) = it 
must hold that I(B(r)) < 1(1). Now, for each / e L and r e R p (L)<~)Pi by definition 
it holds that L n supp(I) = 0, meaning I(B(r)) < 1(1). Thus, sup{/(P(r)) | r € 
P P (L)} < sup{7(Z) I I e L}, meaning I \h h¥(L,P). □ 



Now, we can extend the ASSAT-procedure from (Lin and Zhao 2004) to fuzzy 



answer set programs P. The main idea of this method is to use fuzzy SAT solving 
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techniques to find models of the fuzzy prepositional theory which consists of the 
completion of P, together with the loop formulas of particular maximal loops of P. 
If a model is found which is not an answer set, then we determine a loop that is 
violated by the model and add its loop formula to the fuzzy propositional theory, 
after which the fuzzy SAT solver is invoked again. The algorithm thus becomes: 

1. Initialize Loops = 

2. Generate a model M of comp(I) U LF(P, Loops), where LF(P, Loops) is the 
set of loop formulas of all loops in Loops. 

3. If M = U* pM , return M as it is an answer set. Else, find the loops occurring 
in supp{I Q IIp M ), add their loop formulas to Loops and return to step 2. 

The reason that we can expect this process to be efficient is articulated by Propo- 
sition [6j Indeed, when searching for violated loops, we can restrict our attention to 
subsets of supp(I Q lip,). Although the worst-case complexity of this algorithm is 
still exponential, in most practical applications, we can expect supp(I QH* pI ) to be 
small, as well as the number of iterations of the process that is needed before an 



answer set is found. In (Lin and Zhao 2004) experimental evidence for this claim 



is provided in the case of classical ASP. Last, note that the fuzzy SAT solving 
technique depends on the t-norms used in the program. If only the Lukasiewicz 
t-norm is used, we can use (bounded) mixed integer programming (bMIP) ( |Hahnle| 



1994). Since Fuzzy Description Logic Solvers are based on the same techniques as 
fuzzy SAT solvers, we also know that for the product t-norm we need to resort to 
bounded mixed integer quadratically constrained programming (bMICQP) ( |Bobillo 



and Straccia 2007) 



6 Example: the ATM location selection problem 

In this section we illustrate our algorithm on a FASP program modeling a real-life 
problem. Suppose we are tasked with placing k ATM machines ATM — {oi, • • • , dk} 
on roads connecting n towns Towns = {t\, . . . , t n } such that the distance between 
each town and some ATM machine is minimized, i.e. we aim to find a configu- 
ration in which each town has an ATM that is as close as practically possible. 
To obtain this we optimize the sum of closeness degrees for each town and ATM. 
Note that this problem closely resembles the well-known fc-center selection prob- 



lem (see e.g. (Ausiello et al. 1999)). The difference is that in the fc-center problem 
the ATMs need to be placed in towns, where we allow them to be placed on the 
roads connecting towns. We can model this problem as an undirected weighted 
graph G = (V,E) where V — Towns is the set of vertices and the edge set E 
connects two towns if they are directly connected by a road. Given a distance func- 
tion d : Towns x Towns — > K that models the distance between two townsj^J the 
weight of the edge (a, b) € E is given by the normalized distance d(a,b)/d sum , 
where d sum = Y!,{d(h,te) I h,h e Towns}. 



2 For cities that are not connected the function d models the distance of the shortest path between 
them. 
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Since our FASP programs can only have t-norms in rule bodies, we also need to 
find a way to sum up the distances between towns and ATM machines. By using the 
nearness degree, or closeness degree, which for a normalized distance d is defined as 
1 — d, we can perform summations of distances in our program. To see this, consider 
the following derivation: 

71(1 — disti, 1 — dist 2 ) = max(l — disti + 1 — dist 2 — 1, 0) 
= max(l — (disti + dist2), 0) 
= 1 — min(disti + dist 2 , 1) 

Hence, by applying the Lukasiewicz t-norm on the nearness degrees, we are summing 
the distances. 

The program Patm solving the ATM selection problem is given as follows: 



gloc : 


loc(A,Tl,T2) 


<— 


Ti(conn(Tl,T2),P) 


gnear : 


locNear{A,Tl) 


<— 


M(locNear'(A,Tl)) 


gnear' : 


locNear'(A,Tl) 


<- 


Ti(loc(A, Tl, T2),Af l (near(Tl, T2)), 
locNear(A,T2)),Tl ± T2 


nearr : 


near(Tl,T2) 


«- 


Ti(conn(Tl, T3),near(Tl, T3), near(T'S, T2)) 


locr : 


loc(A,Tl,T2) 




loc(A,T2,Tl) 


atmr : 


ATMNear(A, T) 


<- 


Ti(loc(A, Tl, T2), locNear(A, Tl),near{T, Tl)) 


Wist : 


totNear 


«- 


Ti{{ATMNear(a, t) \ a G ATM , t G Towns}) 


where 









P = Ti({Af m (loc(A,Tl',T2')) | {T[X} {Ti,T 2 }}) 
Note that, due to grounding, a rule such as locr actually corresponds to a set of 
variable-free rules {locr a ^ lt t 2 I a € ATM,t\,t2 € Towns}. We will keep referring 
to the specific grounded instance of a rule by the subscript. 

Program Patm consists of a generate and define part, which for a specific con- 
figuration is augmented with an input part consisting of facts. The generate part 
consists of the three rules gloc, gnear, and gnear' , which generate a specific config- 
uration of ATMs. The gloc rule chooses an edge on which the ATM machine A is 
placed by guessing a location for an ATM that does not yet has an assigned location, 
as ensured by the f3 part of this rule. The gnear and gnear' rules generate a location 
on this edge where A is placed. Rules gnear and gnear' originate from the constraint 
d(a, ti) = d(ti,t 2 ) — d(a, t 2 ), where d(x, y) is the distance between x and y, if ATM 
a is placed on the edge between ti and t 2 . Defining n(x,y) as the nearness degree 
between x and y and noting that n(a,ti) = 1 — d(a,t\) = 1 — (d(ti,t 2 ) — d(a,t 2 )), 
we can rewrite this constraint in terms of t-norms and nearness degrees: 

n(a, t\) = l- (d(h,t 2 ) - d(a, t 2 )) 

= 1 - (d(ti,t 2 ) + (1 - d(a,t 2 )) - 1) 
= l-Ti(d(ti,t 2 ),l-d(a,t 2 )) 
= 1 - 71(1 - n(ti,t 2 ),n(a, t 2 j) 
= M(7T(1 -n(t u t 2 ),n{a,t 2 ))) 
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Hence, the bodies of rules gnear and gnear' ensure that this constraint is satis- 
fied. The reason we need two rules and cannot directly write a rule with body 
Af s (Ti(loc(A,Tl,T2),JVi{near(Tl,T2)),locNear(A,T2)) is that the syntax does 
not allow negation in front of arbitrary expressions. 

Rule nearr recursively defines the degree of closeness between two towns based 
on the known distances for connected towns. Additionally, since the bodies of rules 
with the same head are combined using the maximum, the nearness degree obtained 
by nearr is always one minus the distance of the shortest path. The loer rule 
makes sure that if an ATM is located on the edge between town Tl and T2, it 
is also recognized as being on the edge between T2 and Tl, as we are working 
with an undirected graph. The atmr rule defines the location between a particular 
ATM machine and a town. Note that due to rule loer this rule also covers the 
case when near(T,T2) is higher than near(T,Tl). The tDist rule aggregates the 
total distances such that different answer sets of this program can be compared 
and ordered. In this way we could for example search for the answer set that has a 
maximal total degree of nearness, i.e. in which the distance from the towns to the 
ATMs is lowest. 

Consider the specific configuration Gp = (V,E) of towns Towns — {^1^2,^3} 
depicted in Figure[3]and suppose ATM — {a\, a 2 }. In Figure|4]we depicted a subset 
of the dependency graph of the grounded version of P'atm = Patm U F, where F 
is the input part of the problem, given by the following rules 

F ={conn(t,t') <- 1 | t,t' e Towns, (t,t') G E} 

U{near(t, t') <- k \ t, t' e Towns, (t, t') € E, k = 1 - (d(t, t')/d sum )} 

For the configuration depicted in Figure [3] the input part F is 

F ={conn(ti,ti) <— 1, conn(t\, £2) 1, conn(ti, t^) <— 1} 
U {conn(t2,ti) <— 1, conn(t 2 , £2) 1) connfe, £3) 1} 
U {connfts, ti) <— 1, conn(t%, £2) 1, conn(ts, £3) <— 1} 
U {near(ti, t\) <— 1, near(t\, £2) 4— 0.8, near(t\, £3) <— 0.7} 
U {near(i 2 , t\) <— 0.8, near(t 2 , £2) 1, nearfo, £3) <— 0.5} 
U {near(t3,ti) <— 0.7, near{t^, £2) <— 0.5, nearit^, £3) <— 1} 

It is clear that P'atm contains a number of loops. The completion of P'atm IS the- 
following fuzzy propositional theory: 

conn{t\,ti) w 1, conn{t\,t2) ~ 1, conn(ti,t^) « 1 
conn(t2,ti) R3 1, conn{t2,t2) ~ 1, conn(t2,t 3 ) « 1 
conn(t^,ti) sa 1, conn(t^,t2) ~ 1, conn^^t^) sa 1 
near(ti, t\) « 1, near(t\, i 2 ) sa 0.8, near(t\, t%) sa 0.7 
near(J,2, t\) sa 0.8, near(t,2, ^2) ~ 1, near(t2, £3) w 0.5 
near(ts, t\) sa 0.7, near(ts, £2) ~ 0.5, near{t^, £3) sa 1 
Zoc(ai,*i,*i) max(7I(conn(t 1 ,i 1 ),/3 ljlil ),/oc(a 1 ,i 1 ,< 1 )) 
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loc(a\,t\,t 2 ) ~i max.(Ti(conn(ti,t2), Pi,i,2),loc(ai,t 2 ,ti)) 
loc(a 1 ,t 1 ,t 3 ) max(Ti(conn(t 1 ,t 3 ,/3 1 . lt3 ) 1 loc(a 1 ,t 3 ,t 1 )) 

loc(a 2 ,t 3 ,t 1 ) Rii ma,x(Ti(conn(t 3 ,t 1 ),f3 2 , 3 ,i),loc(a 2 ,t 1 ,t 3 )) 
loc(a 2 ,t 3 ,t 2 ) w ; ma,x(Ti(conn(t 3 ,t 2 ),f3 2 , 3t2 ),loc(a 2 ,t 2 ,t 3 )) 
loc(a 2 ,t 3 ,t 3 ) « ; ma,x(Ti(conn(t 3 ,t 3 ),f3 2 , 3t3 ),loc(a 2 ,t 3 ,t 3 )) 
locN ear{a\,t\) Ni{locNear' (cii, ii)) 

locNear(a 2 ,t 3 ) Ni(locNear' (a 2 ,t 3 )) 

locNear' (ai,ti) max(7I(Zoc(ai, t\, t 2 ), locN ear (ai,t 2 ),Afi (near (t±, t 2 ))), 
Ti(loc(ai,ti,t 3 ),locNear(ai,t 3 ),M(near(ti,t 3 )))) 

locNear' \a 2 , t 3 ) wj max(Ti(loc(a 2 ,t 3 ,ti),locNear(a 2 , ti),Mi(near(t 3 ,t\))), 
Ti(loc(a 2 ,t 3 ,t 2 ),locNear(a 2 , t 2 ),Mi(near(t 3 ,t 2 )))) 

near(t\,ti) ~i ma,x(Ti(conn(ti,ti),near(ti,ti),near(ti,ti)), 
Ti(conn(ti,t 2 ),near(t 1} t 2 ),near(t 2 , ii)), 
Ti(conn(ti,t 3 ),near(ti,t 3 ),near(t 3 , ti)), 1) 

near(ti,t 2 ) max(7z(conn(ii, ii), near(t\, t\), near(ti, t 2 )), 
Ti(conn(ti,t 2 ),near(t 1} t 2 ),near(t 2 , t 2 )), 
Ti(conn(ti,t 3 ),near(t 1} t 3 ),near(t 3 , t 2 )), 0.8) 

near(t 3 ,t 3 ) « ; ma,x(Ti(conn(t 3 ,t 3 ),near(t 3 ,t 3 ),near(t 3 ,t 3 )), 
Ti(conn(t 3 , t 2 ), near (t 3 ,t 2 ), near (t 2 , t 3 )), 
7I(conn(t3,ti),near(i 3 ,ti),near(ti,t 3 )),l) 

ATMNear(ai, ti) s=y max(7z(/oc(ai, ti, ti), locNear(a\, ii), near(ti, tij), 

Ti(loc(ai 7 ti 7 t 2 ),locNear(ai,ti) 7 near(ti,tij), 

Ti{loc(a\,t 3 ,t 2 ) , locNear(a\, t 3 ),near{t\,t 3 )) 
Ti(loc(ai,t 3 ,t 3 ),locNear(ai,t 3 ),near(ti 1 t 3 ))) 

ATMNear(a 2 , t 3 ) w max(7I(/oc(a2, ti, ti), locNear(a 2 , t\), near(t 3 , tij), 

Ti(loc(a 2 ,ti,t 2 ),locNear(a 2 ,ti),near(t 3l ti)), 

Ti{loc(a 2 ,t 3 ,t 2 ),locNear(a 2 ,t 3 ),near(t 3 , t 3 )) 
Ti(loc(a 2 ,t 3 ,t 3 ),locNear(a 2 ,t 3 ),near(t 3 , t 3 ))) 
totNear w Ti{ATMNear(a, t) \ a € ATM, t e Towns} 
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where 

Hi tj ,k = Ti({N' m (loc(a i ,t' j ,t , k )) | {t' v t' k } ± {*,-,**}}) 

Note that e.g. the 1 in the right-hand side of the fuzzy proposition with near{t\, t{) 
on the right-hand side stems from the inputs F we added to Patm- From the com- 
pletion comp(P' ATM ) we can see that an interpretation M satisfying M(near(tl, t2)) = 
1 can be a model of comp(P' ATM ) , which is clearly unwanted as this would overes- 
timate the nearness degrees between towns (i.e. underestimate the distances). For 
example, consider 



M ={loc(a 1 ,t 1 ,t 2 ) 1 ,loc(ai,t 2 ,t 1 ) 1 ,loc(a 2 , t±, i 3 )\ loc(a 2: h^i) 1 , 

locNear(ai, ti) 1 , locNear(a\, t 2 ) , locNear(a 2 , t\f' 75 , locNear(a 2 , £3) 
locNear'{a 2 , ii) a25 , locNear'(a 2 , i 3 ) a25 , near(ti, ti) 1 , near(ti, t 2 )\ 
near(t 2l ii) 1 , near(t±, t 3 ) ' 7 , near(t 3 , ti) ' 7 , near(t 2l t 3 ) ' 5 , near(t 3 , t 2 ) 1 
near (t 2 ,t 2 ) X , near (i 3 , is) 1 , ATMNear{at , ij ) X , ATMNear{a! , t 2 f , 
J 4TMiVear(a 1 , t 3 ) ' 7 , ATMNear(a 2 , ti) ' 75 , ATMNear(a 2 , t 2 ) 0J5 : 
ATMNear(a 2 ,t 3 )°- 75 } 



Note that atoms a for which M(a) = are not included in the set notation, which 
is e.g. the case for totNear. One can easily verify that M is a model of comp(P' ATM ). 
To check whether M is an answer set we compute \ M by repeatedly applying 

\ ATM) 

II(p' -jM , starting from the empty set, until we obtain a fixpoint, and check whether 
M = II/p, w . Performing this procedure, we obtain 



^*( p 'atm) m = oc ( a ii *i! *a) 1 loc(ai,t 2 , ti) 1 ,loc(a 2 ,ti,t 3 ) 1 ,loc(a 2 ,t 3 , ti) , 
locNear(a\, ii) 1 , locN ear(a\, t 2 ) 1 , locNear(a 2 , ii) 0,75 , 
Zoc7Vear(a 2 , ^ 3 ) a75 , ZocAW(a 2 , ti) ' 25 , locNear'(a 2 , i 3 ) ' 25 , 
near(ti, t-y) , near(ti, i 2 ) ' 8 , near(t 2 ,ti) S , near{t\ 1 1 3 ) ' 5 , 
near(t 3 , ti) ' 5 , near(t 2 ,t 3 ) - 7 , near(t 3 , t 2 ) ' 7 , near (t 2 , t 2 Y , near (t^,,^) 1 , 
ATMNear( ai , ) 1 , ATMNear(cn , ATMNear{ai, < 3 )°' 7 , 
ATMNear(a 2 , ii) ' 75 , ATMNear(a 2 , t 2 ) 0J5 , ATMNear(a 2 ,t 3 )°- 75 } 



We can see that IF P , w (neartti, t^Y) = 0.8 ^ M(near(ti,t 2 )), hence M is not 
an answer set of P' A tm- From Proposition^ we then know that there must be a loop 
in suppfMGlITp, w ) = {nearer, i 2 ), near(t 2 , ty)} whose loop formula is violated. 
Looking at the dependency graph, we can see that L = supp(M Q IF P , v M ) = 
{near(ti,t 2 ),near(t 2 ,ty)} contains three loops: L\ — L, L 2 = {near(ty,t 2 )} and 
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Fig. 3. Town configuration for Patm- The weights on the edges denote the nearness 
degrees between towns tt, t 2 and t 3 



Li = {near(t2,ti)}. Their loop formulas are 
LF(Li, P'atm) — 2^(max [nearfa, ij), nearfa, ti)j , max ( Tiiconnfa, t 3 ), 

near fa, t 3 ), nearfa, t 2 )), 0.8, Ti{connfa, t 3 ), nearfa, t 3 ), near(t 3 , t^fj 

lf(l 2 ,f; tm ) = 

Z(max (near fa, t 2 )^ , max (ri(conn(ti,t 3 ), nearfa, t 3 ), nearfa, t 2 )),0.8j 
h¥(L 3 ,P' ATM ) = 

X(max [nearfa, ti)^j , max {l\(connfa, t 3 ), nearfa, t 3 ), near fa, t\)),0.8j 

Clearly, these loop formulas are violated by M, hence following the algorithm in- 
troduced in Section [5j we create a new fuzzy propositional theory comp(P' ATM ) U 
{LF(Xi, P' ATM ),h¥(L 2 , P'atm) '> LF(X 3 , P ATM )}, an( i try to find a model of this new 
theory. Consider then the following model of this new theory: 



M ={loc(ai, ii, t 2 ) , loc(ai, t 2 , t\) , loc(a 2 , t\, t 3 ) , loc(a 2 , t 3 , h) , 
locNear{a 1 ,t 1 ) - 15 JocNear(a 1 ,t 2 ) - 05 ,locNear'(ai,t 1 ) - 85 , 
locNear'(a 1 ,t 2 )°- 95 locNear{a 2 ,t 1 )°- 75 ,locNear(a 2 , t 3 ) - 75 , 
locNear'(a 2 , h) - 25 ,locNear' (a 2 ,t 3 ) - 25 ^earfa,^) 1 ,nearfa,t 2 )°- S , 
nearfa, ii) ' 8 , near fa, t 3 ) ' 7 , nearfa, t\f' 7 , nearfa,t 3 ) ' 5 ,nearfa, t 2 ) 05 , 
near fa, t 2 ) ,near fa, t 3 f , ATMNear{ ai , t t )°- 85 , ATMNear{ ai , t 2 ) - 95 , 
ATMNear( ai ,t 3 ) - 55 , ATMNearfa,^) - 75 , ATMNear(a 2 ,t 2 ) Q - 55 , 
ATMNear(a 2 ,t 3 ) - 75 } 

One can readily verify that this model is an answer set of P'atm ; hence the 
algorithm stops and returns M. 

We could have solved this problem using Mixed Integer Programming (MIP^ 
However, the exact encoding of this problem would be less clear and straightforward 
to write. The reason for this is that in the MIP translation the loop formulas would 



3 Though in general the Godel negation Mm cannot be implemented in MIP, in the ATM example 
we can implement the gloc rules using integer variables. 
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Fig. 4. Dependency graph of Patm 

need to be explicitly represented in the program, while in FASP this is handled 
implicitly. Hence, only the implementer of a FASP system needs to handle these 
loop formulas, not the developer who writes the FASP programs. This is exactly 
the power of FASP: providing an elegant, concise, and clear modelling language for 
representing continuous problems, which, thanks to the results in this paper, can be 
automatically translated to lower-level languages for solving continuous problems, 
such as MIP. 

7 Discussion 

The reader might wonder why we limit our approach to FASP programs with t- 
norms in their body, because at first sight it seems the presented approach is easily 
extendable to arbitrary functions. It turns out that this is not the case, however. 
Consider FASP with the Lukasiewicz t-norm in rule bodies. As mentioned before, 
the completion of such a program, and its loop formulas, are formulas in Lukasiewicz 
logic and are implementable using MIP. Now let us consider FASP where both the 
Lukasiewicz t-norm and the Lukasiewicz t-conorm may occur in rule bodies. At 
first, one would suspect that the loop formulas of such a program would again be 
formulas in Lukasiewicz logic. This turns out to be wrong however. To see this, 
consider the following rules: 

b<-M(a) 
bi-Si{b, b) 
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One can readily verify that in the answer sets of a program containing these rules, 
literal b will be equal to N m {a) (provided that b does not occur in the head of 
any other rule). However, the negation Af m cannot be implemented in MIP, as 
the solution space of a MIP problem is always a topologically closed set (viz. the 
union of a finite number of polyhedra), whereas the solution space of a constraint 
b « Af m (a) cannot be represented as a closed set due to the strict negation in the 
definition of J\f m . This means that as soon as the Lukasiewicz t-conorm is allowed, 
in general, there will not exist a Lukasiewicz logic theory such that the models of 
that theory coincide with the answer sets of a given program. Hence, it is clear that 
the case where other operators than t-norms are used requires a different strategy. 

Finding generalized loop formulas that cover e.g. both the Lukasiewicz t-norm 
and t-conorm is not a trivial problem. To illustrate some of the issues, let us ex- 
amine two intuitive candidates. First, remark that the loop formulas introduced 
in Section [5] eliminate certain answer sets (i.e. they are too strict). Consider the 
following program P: 

a <— Si (a, b) 
b^k 

where k g [0,1]. This program has one loop, viz. {a, b} with corresponding loop 
formula max(a, b) < k. Now note that for k > the value of a in any answer 
set is equal to 1. Hence, the loop formula incorrectly eliminates all answer sets in 
this case. One might think this can be solved by including a condition in the loop 
formula: (max(a,6) < I) V (b > 0). This formula however fails to eliminate models 
that are not answer sets (i.e. it is not strict enough) on the following program: 

a^r m (5 ; (a,6),0.8) 
b <r- Si (a, b) 
b^k 

If k > the unique answer set of this program is {a 8 , b 1 }. However, {a 1 , b 1 } is 
also a model of the completion of this program and satisfies the above loop formula. 

Although again more refined loop formulas can be thought of that handle the 
latter program correctly, we are pessimistic about the possibility of finding loop 
formulas that cover all cases. It appears that such a general solution should be able 
to capture some underlying idea of recursion: one loop may justify the truth value 
of some atom a, up to a certain level, which may then trigger other rules that justify 
the truth value of a, up to some higher level, etc. 

Note that this problem does not occur in classical ASP (or when using the max- 
imum t-conorm), since e.g. a <— b V c is equivalent to a 4— b and a 4— c, which is 
indeed why disjunctions in the body of rules are not considered in classical ASP. 
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8 Related Work 



The approach to fuzzy answer set programming for which we provided the trans- 
lation to fuzzy SAT is called an unweighted implication-based approach. There also 



exist weighted implication-based approaches (e.g. (Lukasiewicz and Straccia 2007a 



Madrid and Ojeda-Aciego 2011||Madrid and Ojeda-Aciego 2 008 Ma drid and Ojeda- 



Aciego 2009)), which use rules of the form 

r :a^T(b 1 ,...,b n ) (11) 

where r is a rule label, a is an atom, for 1 < i < n, are extended literals, and 
a G [0, 1]. An interpretation / models this rule iff 

l r (T(I(h),...,I(b n )),I(a))>a 

Since T r is the residual implicator of T this is equivalent to 

1(a) >T(/(&i),... ,/(&„), a) 

Hence a weighted rule of the form ( |TT| ) above can be simulated by the rule 

r' :a<r- T(h, . . . ,b n ,a) 

Thus, the translation presented in this paper can equally be applied to weighted 
implication-based approaches. 

In addition to the implication-based approaches (IB) one also finds annotation- 



based (AB) approaches (see e.g. (Straccia 2006)). In the annotation-based setting 
a rule is of the form 



A :/(/?!,...,/?„) <-B 1 



i B n . 



Such a rule asserts that the value of atom A is at least /(/3i, . . . , (3 n ) if the value 
of each atom 1 < n, is at least pV In this setting / is a computable function 
and Pi is either a constant or a variable ranging over an appropriate truth domain. 
Due to the difference in semantics between the IB and AB approaches, our method 
is not directly applicable to AB frameworks. One can find an in-depth overview of 



logic programming with fuzzy logic in (Straccia 2008). 



In (Van Nieuwenborgh et al. 2007a), an implementation method for FASP pro- 



grams with a finite truth value set is presented, which consists of translating a 
FASP program to a specific DLVHEX program. For solving continuous problems, 
however, we need infinite truth values, for which a solving method is much harder 
to construct. Our method is able to handle continuous problems, and additionally 



is more flexible than (Van Nieuwenborgh et al. 2007a I since any method for solving 



continuous problems can be used as the backend, including fuzzy SAT solvers and 
the vast body of existing MIP solvers. 

Apart from fuzzy answer set programming, in recent years possibilistic and prob- 
abilistic answer set programming have been developed. Both of these approaches 
can be reduced to classical SAT. In the case of probabilistic ASP, there is a direct 



translation method (Saad 2009b), while a possibilistic ASP program can be trans- 



lated to an equivalent ASP program, on which the ASSAT procedure can then be 
applied. 
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9 Conclusion 

In this paper we have focused on the translation of FASP programs to particular 
satisfiability problems. We have introduced the completion of a program and have 
shown that in the case of programs without loops, the models of the completion 
are exactly the answer sets. Furthermore, to solve the general problem, we have 
generalized the notion of loop formulas. This translation is important because it 
allows to solve FASP programs using fuzzy SAT solvers. Under appropriate re- 
strictions, for example, the satisfiability problems that are obtained can be solved 
using off-the-shelf mixed integer programming methods. From an application point 
of view, this allows us to encode continuous optimization problems in a declarative 
style which is similar to traditional answer set programming. This style of encoding 
problems is often more intuitive, as well as more concise, while the results we have 
presented ensure that the power of mathematical programming techniques can still 
be employed to find the solutions. 
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